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Abstract 

We revisit the QCD+QED axial anomaly low-energy theorems which give an exact relation 
between the matrix elements of the gluon and photon parts of the axial anomaly operator equation 
within the framework of the effective action derived from the instanton vacuum. The matrix 
elements between the vacuum and two photon states and between the vacuum and two gluon 
states are investigated for arbitrary Nj in the chiral limit. Having gauged the effective action 
properly, we show that the model does exactly satisfy the low-energy theorems. 
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1. The quantum fluctuations may destroy the symmetries of the initial classical La- 
grangian. One of the most important examples is the famous axial anomaly in gauge the- 
ories. The axial anomaly leads to many interesting nonperturbative features intimately 
related to the topologically nontrivial structure of the vacuum. 

The axial anomaly in the divergence of the singlet axial-vector current in QCD + QED 
brings forth a low-energy theorem for the matrix elements of this operator equation between 
the vacuum and two-photon states (see, for example, a review P): 







q 2 ~ 
N f ^—GG 



at q 2 <C m 2 , [m^ is the r/-meson mass around 1 GeV). Nf is the number of flavors, g is the 
QCD coupling constant with GG = \^ vXa G a ^ v G\ a and N c is the number of colors. Here, 
G a ^ y stands for the operator of the gluon field strengths. Cf represent the electric charges of 
quarks. The photon part of Eq.((TJ) can be expressed explicitly as = e^q^—e^q^, where 
e^ 1,2 ) and q\ :2 are, respectively, polarizations and momenta of photons, and q = q% + q 2 . In 
perturbation theory, it leads to at least results of order e 2 g 4 for the left-hand side of Eq. (JTJ). 

Another nontrivial low-energy theorem concerns the matrix element between the vacuum 
and two-gluon states: 

(0|/GG|2gluons) = (2) 

at the same limit as in Eq.(fTj). Hence, the solution of these theorems given by Eqs. (jll2|) 
is only pertinent to the nonperturbative phenomena related to the topologically nontrivial 
structure of the QCD vacuum. 

Without any doubt instantons represent a very important topologically nontrivial com- 
ponent of the QCD vacuum. The most important parameters of the QCD instanton vacuum 
are the average instanton size p and inter-instanton distance R. Shuryak 0, HI estimated 
these two parameters phenomenologically as p ~ 1/3 fm, R ~ 1 fm. It was confirmed 
by theoretical variational calculations jUll and recent lattice simulations of the QCD vac- 
uum 0]. In particular, the spontaneous breakdown of chiral symmetry is realized very 
well via the instanton liquid model. Hence, instantons play a pivotal and more significant 
role in describing the lightest hadrons and their interactions, compared with confinement 
forces 

Some years ago, one of the authors showed that the low-energy QCD effective chiral 
action from the instanton vacuum j^, 0] satisfies the low energy theorems with an accuracy 
of about 17 % 11 1 . This discrepancy is due to the fact that the vector and axial- vector 



currents are not conserved because of the momentum-dependent dynamical quark mass in 
the effective action [l2l . fl3l 0|. The nonconservation of the Nother currents breaks 
various Ward identities in the model. In order to remedy these problems, we shall gauge 
the low-energy QCD effective partition function properly, from which the conserved currents 
can be derived. 

In the present paper, we shall continue to test the low-energy QCD effective chiral action 
from the instanton vacuum in the chiral limit aiming at resolving the discrepancy existing 
in Ref. ^(J. Having gauged the low-energy QCD effective partition function, we shall show 
that it does exactly satisfy these theorems as it should be. 

2. We have to calculate the quark determinant in the presence of instantons and external 
electromagnetic field v M , and then average it over the collective coordinates of the instantons 



2 



to find the effective low-energy _QCD partition function. The first thing we have to do is to 
find the extended zero-modes $o for the quark placed into the instanton field A 1 ^ and the 
electromagnetic field v^, which obey the following equation: 

(ip x + gA\x)+ei,(x))%{x) = (3) 
The extended zero-modes satisfy the U(l) gauge transformation 

$o = -^0) L(x) = exp ( iex^ / v IM (ax)da ) , (4) 



where L is simply a path-ordered exponent Pexp(i Jq v u (s)ds u ) with a straight path con- 
necting the points to x [3, 17 1. Then Eq. © is reduced to 

{ip x +g/i 1 (x)+ef(x))& (x) = 0, (5) 
where 

v 'ni x ) = x p I F Pfl (ax)ada (6) 

J 

has an explicit gauge-invariant form. An approximate solution of Eq. (JHJ is 

$' = $ - S'ef^o, (7) 

which can be expanded to any desired order of ev'. Here, S 1 = (ifi + gA 1 + im)~ l is a quark 
propagator (with a small mass m) in the field of a single instanton: 

j _ goKgo] _ ^ |$n)($n| ra v 

<-> — : r onz, &NZ — y, -. — — \p) 

im ^ \ n + vm 

where $ n is defined by {%</) + gA 1 )^^ = A n |$ n ) (with the convention A = 0). The quark 
propagator of the non-zero mode Sjsrz in the single instanton field was derived in an exactly 
closed form [18]. It is clear from this formula that the quark propagator is reduced to the 
free one at a short distance as well as a long distance. We follow here the argument by 
Ref. |H and will assume the following relation: 

S NZ {x,y)~S {x-y). (9) 

Combining Eq.Q with Eqs. (j4l5l7|) . we obtain the extended zero mode $o = L& . Now, we 
consider the quark propagator S in the presence of the single instanton A 1 ^ and with the 
assumption given by Eq.©, so that we obtain for the quark propagator S 1 : 

S 1 = (ip + gA 1 + ef + im)- 1 = S 1 - S^p 1 + ■■■ 

- !*•> <*»l +Snz -( JMM + Snz ] ev ( + Skm ) + ... 



im \ im / V im 



& + (10) 



im 



where So = (ifi + ef> + im) -1 . 

Now, we are in a position to consider the ensemble of instantons A^ = J2i A?>> represent- 
ing the instanton-liquid picture of the QCD vacuum. This sum consists of both instantons 
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and antiinstantons. The quarks are moving in the presence of the instanton ensemble and 
external electromagnetic field v. Thus, we can derive an extended zero-mode approximation 
for the quark determinant. 

The quark propagator S can be expanded with respect to a single instanton: 

S = S - J2 SogAiSo + E SogAjSogAjSo + . . . (11) 
/ i, j 

= So + - So) + £(& - S^So^Sj - So) (12) 

+ £ (Si — Sq)Sq (Sj — So)S (Sk — Sq) + . . . . 

Here, the indices I,J,K,--- designate both instantons and anti-instantons. In Eq.(|12|). 
we have resummed the contributions corresponding to an individual instanton. Sj stands 
for the quark propagator in the fields of the single instanton I as well as in the external 
electromagnetic field v and represents the sum of the zero and non-zero modes. At a large 
distance from the center of the instanton, Si becomes the propagator So in the presence of 
the external field v. Equation ()12|) is obviously an extended quark propagator of that from 
Ref. (see also review 1^|). Thus, Eq. (JTUJ) can be written as 

(5/ - So) (ar, y) ~ : ! . (13) 

v ' im 

In the present case, the expansion given by Eq.([12p becomes 



! im 



+ V^^f [d 4 r^ o (r)(^ft + e0 + lm)$J fi (r))^^ + ■■■. 
jj-j im \J ' J im 

Defining the following overlapping integrals 

o.u = - j d*r^ i0 (r)(^ + ^)$ J>0 (r), (15) 
we are led to the expression: 

S(x, y) ~ S (x, y) + £ <f> 7 , (x) ( ) & J>0 (y). (16) 

\a + im/ij 

The total quark determinant should be now splitted into two parts, i.e. low and high 
frequencies (with respect to a mass parameter Mi ~ p _1 ): Det = Deti ow x Dethigh 0]- The 
high-frequency part Dethigh can be written as a product of the determinants in the field of 
individual instantons, while the low-frequency one Deti ow has to be treated approximately, 
would-be zero modes being taken into account only. Thus, we get for the low-frequency part 
of the quark determinant as follows: 

lnDet low = ]T / idm[Tr(S - S ) - Tr(S - S )\, (17) 

J JMi 
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where 

Tr(S - So) = Tr^— . (18) 
a + im 

The second term in Eq. (|17jl can be included in the normalization factor. Taking into account 
the external gauge field v, we get 

Deti ow = det^v = det B, Bjj = djj + im5jj. (19) 

If we had turned off the external field v, we would have obtained exactly the same quark 
determinant det at by Lee and Bardeen [20]. Having fermionized the representation of 



Eq. (|T9*|) |l0i lllj , we derive the gauged quark determinant in terms of the constituent quarks 
det N , = J DipD^ exp( J d 4 x ^ ip\(ifi + erftyf) 

N+ N- 

x u(m im f-vM^f])U( im f-v-ii>h^f})i (20) 
/ + 



where 

V±[1>\,il>f] = I d 4 x(i>}(x)(ip + erf)$±,o(x;Z ± )) J A(^L,o(y; t±)W + erfMv)) (21) 
d 4 x(^ f (x)L f (x)tp^ ±fi (x; e±)) / d*y{& ±t0 (y; ^)ipL){y)^ f {y)). (22) 



Here, we have used Eqs. (l7J9j) to derive Eq.(|22j) from Eq.(|21|l. The range of the integration 
in Eq.(j2*Tj) is cut at p, since it is defined by zero-mode functions $±,o- In the present case, 
we assume the soft-photon external field qp <C 1 (q is a photon momentum) which simplifies 
the gauge factor: 

L f (x) = exp (ie f (x - z) M v^{z + a(x - z))do^j ~ exp (ie f (x - z)^^)) (23) 

Note that external t> M field gauges not only the kinetic term of the effective action but 
also its interaction term VifV'/, tpf] in Eq. (|21|). The reason is obvious: It is the nonlocal 
interaction induced by instantons. The external t^-field is present here due to the factor Lf 
attached to each fermionic line. This factor provides us an approximate gauge invariance 
of the interaction term V±[ipf,ipf] under the gauge transformation which is valid in the 
soft-photon limit qp <C 1. The remaining problem is to average the quark determinant 
over collective coordinates £±. It is a rather simple procedure, since the low density of the 

instanton medium (tt 2 (j^) ~ 0.1) allows us to average over positions and orientations of the 
instantons independently. In the following, we will assume the chiral limit. Thus, Eq. (j20j) 
gives us the gauged partition function: 

Z N [v] = J D^D^exp ^ J d 4 x X>/0'^ + e /^/j ^+ + ^-~> ( 24 ) 
where 

r ~ ( 4-Tr 2 n 2 \ Nf r d 4 ? 



/ 
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and 

d 4 kd 4 l 



J+iz. 



fa 



PkdV i i 

— exp(-*(*-Qs) F{{k+e f v{z)f )F{{l+e g v^)) 2 ) ^)-(l± 7 5 W0-(26) 



The form factor F(k 2 ) is related to the zero-mode wave function in momentum space 
$±(&;£±) and turns out to be F{k 2 ) = -tf t [I (t)K (t) - h^K^t)} , t = \\fkp-p, where 
Ii and Ki are modified Bessel functions. Thus, we find that the gauged partition function 
in the soft-photon limit is nothing but the original partition function with the gauged form 
factor F((k + efv(z)) 2 ). Details of how to derive the effective action from the partition func- 
tion are well described in Ref. 0. The self-consistency condition of the exponentiation and 
bosonization at the common saddle-point leads to the following equation for the dynamical 
quark mass M: 

JAr f d 4 k M 2 F 4 (k) N 

AM / i_l = 27) 

C J (2vr) 4 M 2 F 4 (A;) + k 2 V' 1 ' 

where N/V denotes the finite density of instantons and antiinstantons in Euclidean space. 

In the quasiclassical (saddle-point) approximation, any gluon operator receives its main 
contribution from the instanton background. As an example, for one instanton (anti- 
instanton)J(J), the following gluonic operators can be expressed in terms of the instanton 
background: 

9 2 G 2 {x)= - 192p \ =f( x -z), g 2 GG{x)=±f{x-z). (28) 



[p 2 + (x — z] 



The derivation of the correlators with the gluonic operator GG comes down to that of the 
following partition function: 

Z n [v,k] = Z N l J D^D^exp(-S eS ), (29) 

where the effective action, S e g in the presence of an external fields v^x) and k(x) coupled 
to GG, is given by 



.S.,f = / d 4 x ^fiip + ef^f+i^—j J d 4 z det(iM + (z)J+(z)) 



/2V\ N f~ 1 r 

+ f — J J d 4 z det(zAf_(z)J_(z)), (30) 

where M±(z) = (1 ± J dxK(x)f(x - z)) (Ar/ ~ 1)_1 M 0. Using the bozonization again, one 
can show that the effective action S e s[Ai±, v, k] turns out to be 

2V" 



S cS [M±,v,k} = Jd 4 zl-(N f -l)^—^j (detMi)^- 1 



Tr In 



P + iMF{P 2 )M + F{P 2 ) (1 + («/)) V i(l + 75 ) 



+ iMF{P 2 )M-F{P 2 ) (1 - («/)) V 1(1 - 75 ) 



(31) 
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in the presence of the external fields and k. Here, P M denotes the covariant momentum 
Pfi '■= Vti + zQvn with eQfg = Cf5f g . The matrices M.± describe the meson fields and 
(«/) = I d 4 xK(x)f(x - z). 

Since the matrix element between the vacuum and two-photon state is defined by the 

we need to keep only the terms of order 0(k) 



functional derivative 



<5Zjv [k,v] 



k,v=0 



and 0{v 2 ) in S eS [M± = l,v, k] defined by Eq.fjSU). 

We will employ for the explicit calculation the Schwinger method 2l|, which was devel- 
oped for quantum electrodynamics and later extended to QCD 22l . in order to expand S c s 
with respect to derivatives of the photon background field. Setting A4± = 1, considering 
the first order of k, and defining M(P) = MF 2 (P), we get 



S eS [M± = 1,v,k] = Trln[P + iM(P)\ + — Tr \(P + iM(P))- 1 M(P)(Kf)-y 5 

jyif 



(32) 



where Tr = tr c / 7 J d 4 x(x\ \x) and tr c / 7 denotes the trace over color, flavor, and Dirac 

spaces. The quark propagator with the covariant derivative can be treated as follows: 



(P + tM(P))- 1 = ( P 2 + M 2 (P) + e -Q-o ■ F + M(P)]\ 1 (P - iM(P)), 



(33) 



where a ■ F = a^F^. We need at least four 7 matrices with one 75 so that we may obtain a 
nonvanishing result with the trace over Dirac space (tr 7 ). Thus, the relevant nonvanishing 
part of the effective action is: 



[P 2 + M 2 (P))- X [^-a ■ F + i[p, M{P)\ 



x (P 2 + M 2 (P))-\M 2 (P) + zPM(P))( K f) l5 } . 
We observe that [P, M(P)] is of order 0(e) and M is a function of P 2 . Then 



(34) 



S k f = — Tr 



N f 



r2 (p 2 ))- lfeQ 

-i e Q 



(p 2 + M (p )) [-ftr-F) ( P 2 + M 2 (p 2 ))-'M 2 (p 2 )(Kf) l5 



\p 2 + M\p 2 )Y^a • P, (p 2 + M 2 (p 2 ))^[P, M(P 2 )]} 
x (p 2 + M 2 (p 2 ))-^M(p 2 )(kf) l5 \ +0(e 3 ). 



(35) 



Here, we have taken into account the term to order 0(e 2 ) only. Hence, what we need is to 
calculate [P, M(P 2 )], which becomes in the gauge [p M ,f M ] = 0: 



[P, M(P 2 )] = -2ieQF^ lvPtl 



dM(p 2 
dp 2 



0(e 2 ,d 2 ). 



(36) 



Putting them together, we finally arrive at the expression: 



S, 



hf 



d±p M 2 (p 2 ' 



V dp 2 



(2tt) 4 (p 2 + M 2 (p 2 )f 



tr(^x-P) (kf) l5 . 



(37) 



It was shown in the previous paper [T^j that if we neglect the momentum dependence of 
the M, then we exactly reproduce the low-energy theorem in Eq.(JJ). Hence, In order to 
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prove that the gauged effective action exactly satisfies the low-energy theorem of the axial 
anomaly, it is enough to show that the following ratio R = with 

becomes unity. Replacing the variable |l4| by s = M - , we immediately have the ratio 
R = 1. As a result, the nonlocal chiral quark model from the instanton vacuum exactly 
satisfies the low-energy theorem given by Eq.(jT|). 

Now, we briefly present some calculations related to Eq.(j2J). This matrix element can be 
written in the form: 

(0\g 2 GG\g^ 1 \q 1 ),g(^,q 2 )) 



,(l)oi (2)aa / f)2f)2 

Ml C M2 



| d\d\ (0\Tg 2 GGA';\{x 1 )A^{x 2 )\0)expt{q 1 x 1 + q 2 x 2 )dx 1 dx 2 , (39) 



where denotes a total gluon field. e^ ai and $ are the polarization vectors and the 

momentum of gluons, respectively. 

As usual, we expand the total field (x) around the instanton background. The main 
term in Eq.()39|) is the contribution of the instanton background and is of order (D(g~ 2 ). 
The next term is due to the perturbative fluctuations over the instanton background and 
corresponds to the contribution of order 0(g 2 ). It is easy to see from previous considerations 
that the term to order 0(g~ 2 ) is given by the following expression: 
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N 



J ItyDtf exp(-S cfr ) (Y GGAA+ (x) + Y GdAA _{x)) (40) 



with 

;2V\ N f- 1 

Yggaa± = ^ 



(— ) ~ {iM) N * Jd'zf(x-z) 
x / dO(-aD4f )ai (^i)(-9 2 2 )4( J > 2 (^)det J ± (^), (41) 

where the instanton(anti-instanton) is located at the point z with its orientation O 0] over 
which we integrate. 

Repeating the bosonization, we obtain the result for the 0(g~ 2 ) contribution which is 
proportional to TY[(iP + iMF 2 )~ l iMF 2/ y 5 ]. Hence, one can easily show that the contribution 
from the term to order (D(g~ 2 ) vanishes. 

The contribution from the next order (0(g 2 )) comes from two different diagrams. The 
first diagram is the direct contribution of the operator g 2 GG which is equal to —g 2 G t - 1 '& 2 \ 
where 2GWG( 2 ) = e» uX<T G$ a G ( £ a , Gf u a = efq iv - efq^. The factors at the vertices of the 
second one are g7 M A a /2 and iMfF 2 'j 5 Ny 1 . 

A comparison with previous calculations ends up with the result that the contribution 
from the second loop-diagram is equal in magnitude but opposite in sign to that from the 
first one at q 2 — > 0. Because of this cancellation, the contribution from order 0(g 2 ) vanishes 
in the limit q 2 — » 0. Hence, we arrive exactly at the same conclusion as in the previous case. 

From the above calculations, we conclude that the nonlocal chiral quark model from 
the instanton vacuum does exactly satisfy the low-energy theorems, once we have properly 
gauged the low-energy QCD partition function from the instanton vacuum. 
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3. The solution of the low-energy theorems of the axial anomaly is nontrivially related to 
the nonperturbative instanton structure of the QCD vacuum. The effective action approach 
based on the instanton vacuum exactly satisfies these low-energy theorems. 

This approach provides a solid ground for future investigations into different amplitudes 
for the nonperturbative conversion of gluons into hadrons and photons. Further studies 
relevant to all of these problems are under way. 
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